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ABSTRACT 

This paper describes a technique of obtaining particle size 
information from holograms of combustion products. The holograms 
are obtained with a pulsed ruby laser through windows in a com- 
bustion chamber* The reconstruction is done with a krypton laser 
with the real image being viewed through a microscope* The parti- 
cle size information is measured with a Quantimet 720 image 
processing system which can discriminate various features and 
perform measurements of the portions of interest in the image* 
Various problems that arise in the technique are discussed, 
especially those that are a consequence of the speckle due to the 
diffuse illumination used in the recording process* 


INTRODUCTION 

A continuing investigation^] is being conducted at the 
Naval Postgraduate School (NPS) to obtain quantitative data that 
can be used to relate solid rocket propellant composition and 
operating environment to the behavior of solid particulates (Al, 
AI 2 O 3 ) within the grain port and exhaust nozzle* These data are 
needed to improve predictive capabilities of propellant perfor- 
mance, to provide input related to ammonium perchlorate (AP)— 
aluminum interactions for steady— state combustion models, and to 
provide in— motor particle size distributions to allow more accu- 
rate predictions of damping in stability analyses* The expected 
range of particle sizes is from 1 |im to 200 pm* While work has 
been done on a variety of techniques, including high speed motion 
pictures, analysis of post-fire residue with a scanning electron 
microscope, and particle sizing by scattered light measurements, 
this paper will focus on the application of particle sizing using 
holograms of burning propellant strands and slabs in a cross— flow 
environment* 


HOLOGRAM RECORDING 

A diagram of the two-dimensional motor is shown in figure 1. 
A glass window (18.5 mm diameter) in the motor walls allows 
recording of the hologram. The window position allows centering 


♦This research was sponsored by the Air Force Rocket Propulsion 
Laboratory, Edwards Air Force Base. 


PRECEDING PAGE BLANK NOT FILMED 


589 


mm 


of the edge of the largest combustion slab. The smallest combus- 
tion slabs are out of view with these window positions. Self- 
pressurization of the burning slabs was inadequate to achieve 
desired values (34 atm. and greater), so a nitrogen pressuriza- 
tion technique was incorporated into the top of the motor cham- 
ber. A 0-80 screw was mounted outside of the rear window (from 
the laser) to provide a size and position reference in the recon- 
structed hologram. Additionally, holograms were recorded of reso- 
lution charts and glass beads (both opaque and transparent) of 
known size for calibration purposes and investigations in chang- 
ing recording or reconstruction geometries. 

The laser system is a pulsed ruby laser[2] consisting of a 
Q- switched oscillator, a ruby amplifier, beam expanding optics, 
an alignment autoeollisator, a EeNe pointing laser, the capacitor 
bank, and associated electronic power supplies. The operating 
wavelength is 694.3 nm with a beam diameter of approximately 3.2 
cm. A one joule pulse with a 30 ns pulselength was used for this 
investigation. 

All holograms were recorded with diffuse illumination from 
the laser in order to minimize the presence of Schlieren inter- 
ference fringes produced by temperature and density variations of 
the combustion gas products during the burn. The diffuse illumi- 
nation was made by introducing a glass diffuser into the scene 
illumination beam. This diffuse illumination introduces speckle 
into the reconstructed images as will be discussed in detail 
later. The primary problem in sizing the particles in the recon- 
structed image is interference from this speckle which can have a 
maximum size that is comparable to the particles at the lower end 
of the expected particle size distribution. 

The holocamera [3] (shown in the reconstruction diagram of 
figure 2 ) incorporates assisting lenses to increase the effective 
field of view of the camera with a subsequent increase in resolu- 
tion. Holographic plates from AGFA-Gevaert (Model 8E75 HD plates) 
were used for hologram recording. 

Holographic recordings have been made successfully using 
propellant strands burned at operating pressures of 34 and 68 
atm. and with various concentrations of alumininum (up to 15% 
aluminum). Single pulse holography provides a means of effective- 
ly stopping motion. It only provides information during the 
single instant of time, however. Smoke generation (i.e., small 
AI 2 O 3 and binder products, etc.) presents a major obstacle to 
obtaining good holograms. The smoke reduces illumination levels 
in portions of the image and also adversely affects the scene-to- 
reference illumination levels on the hologram recording medium. 
Minimization of the effects of smoke requires experimental de- 
termination of the optimum propellant geometry and the optimum 
delay time for recording the hologram during the burn. 
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HOLOGRAM RECONSTRUCTION 


During reconstruction* the developed hologram plate is re- 
attached to the plate holder and placed into the holocamera box 
(which is mounted on a digitally controlled xyz stage). The plate 
is illuminated with a krypton-argon laser at an angle of approxi- 
mately 60° from normal* as shown in figure 2. The krypton laser 
operates at 647.1 nm, causing slight lateral demagnification of 
the image equal to the ratio of the reconstructing and recording 
wavelengths. A longitudinal distortion equal to the square of the 
wavelength ratio is also expected. The reconstruction laser wave- 
length was selected close to the ruby laser wavelength to mini- 
mize these effects. The real image of the hologram was focused 
onto a spinning mylar disk that was introduced to reduce speckle 
effects in the observed image. The spinning disk changes the 
speckle pattern at a rate faster than the response time of eye or 
imaging system causing a reduction in the contrast of the speckle 
pattern. A variable power microscope was used to view the recon- 
structed image either by eye or with the image scanner of the 
image processing system. Reconstructed images can be recorded by 
a 35mm camera from either the screen of the image processing 
system or from a camera attachment on the microscope. Photographs 
of typical reconstructed images are shown in figure 3. 


IMAGE PROCESSING EQUIPMENT 

The Quantiment 720[4] is a general purpose television-type 
image analyzer that is capable of elementary shape recognition 
and various physical measurements of objects by distinguishing 
differences in grey levels in the image and performing various 
logical tests on measured dimensions. The system is designed in a 
modular fashion so that additional capabilities can be acquired 
with extra electronic modules. When used or configured as a 
'basic' Quantimet 720, the operator must provide all control and 
direct all operations. In the 'advanced mode', a computer is used 
as a controller where all operations are programmed from prior 
experience. The Quantimet at NPS contains the necessary modules 
to be operated in the advanced mode but is currently being uti- 
lized in the basic mode. 

Figure 4 is a block diagram of the image processor in the 
enhanced basic mode. All switch settings must be made on an 
interactive basis and data readings are made by the operator from 
the digital display on the screen. 

The Scanner sends a video signal of of the image to the 
System Control Module where it is digitized and sent to the 
Variable Frame and Scale module. This module* when activated* 
enables an operator to select only a portion of the image for 
analysis, disregarding the rest of the image. The portion of the 
image of interest to the operator is forwarded to the 1— D Auto 
Detector module. This 1-D Auto Detector module can differentiate 
data based on grey level. There are 64 grey levels in an image 
and differentiation can be made of all levels above a certain 
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level, below a certain level, or between two levels. The thresh- 
old levels are under the operator control. Only data that passes 
the desired grey level test is passed to the other modules for 
further processing. The Light Pen module is used by the operator 
to select individual features for further measurements. The Frame 
Smasher module is of use only for long thin objects and is usual- 
ly bypassed in our operations. The Standard Analyzer module can 
measure the selected features of the image. It can measure area, 
perimeter, largest horizontal dimension, largest vertical dimen- 
sion, and other quantities. (All measurements are in units of 
pixels and image calibration is required to reduce the measure- 
ments to physical units.) The Standard Analyzer module can also 
be used to select features for further analysis based on the 
results of these measurements. For example, the analyzer can be 
used to select ail image features whose horizontal (or vertical) 
size exceeds a certain value, whose area exceeds a given value, 
or whose dimensions (either vertical or horizontal) are less than 
a given value. The logical tests based on measured size can be 
combined with the logical tests performed by the 1-D Auto Detect- 
or module for further feature discrimination. The features that 
meet the combined specifications are sent to the Function Compu- 
ter module and the Classifier/ Collector module. These modules can 
perform more measurements on the features of interest or can 
perform still more logical tests on the dimensions. In the cur- 
rent configuration these units are not brought to bear on the 
feature measurement or discrimination problem. 

In the 'automatic' mode (Fig. 5) a computer controls the 
setting of the module functions and the various threshold levels 
of the logical tests, as well as the recording of the data. 
Individual image frames can be transferred into the computer for 
data recording or subsequent analysis. 

A standard image frame is 688 by 880 pixels. As described 
above, a subportion of this frame can be analyzed. The Scanner 
uses a Plumbicon tube selected for its uniformity of response 
over the complete face of the tube. The usual scan rate is 10 
s . Scan rates of 1 s - * and 0.4 s - ^ are switch selectable. The 
display screen is a long-persistence yellow phosphor that can 
accommodate the 10 s ^ scan rate with a minimum of image loss 
before refreshing. The screen can present the raw image, the 
image as processed by the 1— D Auto detector, the image as proces- 
sed by the Standard Analyzer, or any combination of these images. 
A digital readout on the top of the screen presents the quantita- 
tive data to the operator. 


PROBLEM AREAS 

In this section we will describe several identified problems 
in this application of data reduction and potential solutions to 
them. The problems can be broken up into two parts: those that 
exist without the speckle being present and those that exist with 
speckle being present. 


592 


The first problem is to compute the number of locations that 
must be investigated in the hologram (i.e., we must compute the 
number of planes of interest that exist in the volume of the 
reconstruction). The technique of particle sizing requires that 
the system locate the planes of focus of the particles so that 
the particles can be accurately sized. One of the calculations of 
interest is to determine how many planes of best focus will exist 
in the hologram reconstruction. We will assume a nominal depth Az 
of 2 mm to the hologram reconstruction. The operating wavelength 
X is that of the krypton laser — 647.1 nm. A particle of diameter, 
d, will have a depth of focus, Ap, given by 

d 2 

Ap- ( 1 ) 

X 

where the depth of focus is the distance over which the image of 
the particle will be 'in focus'. (This criterion is that the 
diameter of the image of the particle be within a factor of 
(2)^^ of its minimum in-focus value). The minimum number of 
planes, N, that must be investigated to count and size particles 
is given by 


Az 

N » (2) 

Ap 


Az X 



Table 1 gives the value N for particles of differing diameters 
within the range of particle diameters expected in the experi- 
ment. As seen from the table the computational burden is not too 
great for particle sizes larger than 10 pm. As one goes below a 
size of 10 pm, however, the number of sample planes required 
increases rapidly. This is the region where speckle size also 
becomes important, as the size of the speckle is on the same 
order of magnitude as the diameter of the particle. 

Stanton, Caulfield, and Stewart [5] present a technique for 
reducing the computational burden of finding the in— focus 
particle and increasing the accuracy of particle location by 
measuring the Fresnel ring structure of an out-of-focus spherical 
particle. This technique is not available to us since the Fresnel 
rings of the particle can not be seen against the speckle back- 
ground. The particles are not identifiable in the image until the 
observation plane is close to the particle location. 

The Quantimet 720 presents limitations which can lead to 
inaccurate results. One technique of separating particles from 
background is to accept for measurement only those pixels that 
that are darker than a certain threshold. Once the particle have 
been isolated, the particle size is measured. The 720 locates an 


593 


4 H + 


d 

1 Ap 

i 

1 N 

1 

1 

1 

1 

t 

1 pm 

1 1.54 pm 

1 1,294 
1 

““T 

1 

1 

1 

M 

0 

1 

1 154 pm 

1 12.9 

1 

1 

1 

1 

o 

1 3.86 mm 

1 

1 1 
1 

1 

1 

100 pm 

1 1.54 cm 

l 

1 1 

1 

1 



1 

9 

1 

B 

200 pm 

j 6.18 cm 

1 l 

1 


H -I H + 


Table 1 Depth of focus, Ap, and number of independent 
viewing planes, N, in a 2 mm deep hologram volume 


edge of a particle by noting regions where there is a sharp 
change in grey level. The edge location is determined by finding 
the pixels with greatest and least grey level and locating the 
edge at the pixel midway between these two locations. This tech- 
nique is consistent only for particles with the same extremes in 
grey levels. For frames that include particles of differing grey 
levels, the edge location is variable. Uneven illumination is the 
primary cause of differences in particle grey level. The illumin- 
ation problems can occur in the hologram recording (e.g., laser 
nonuniform beam pattern, smoke), in the development of the holo- 
gram leading to a plate with variable transmissivity (bleaching 
the plate would help this problem), and in the reconstruction 
process (e.g., uneven laser illumination, nonuniform spatial 
response of the image tube). Additionally electronic noise in the 
various modules and quantization noise can corrupt the image 
signal leading to measurement errors. 

For the case of having speckle present as in our hologram 
reconstructions (due to the diffuse illumination required to 
eliminate the Schlieren fringes in the reconstruction as previ- 
ously described), the initial point is to reduce the maximum size 
of the expected speckle to the smallest value possible. Speckle 
is the random interference pattern from a diffuse source inter- 
fering with a reference wave. The random pattern of light will 
contain all grey levels from maximum brightness to total black. 
The speckle can have two effects on the image. It produces black 
spots in the background which cannot be readily distinguished 
from real particles. The second effect is that the speckle can 
give the perimeter of the particle a 'swiss cheese 1 appearance 
where the speckle overlaps the edge. This alters any calculated 
measurements of the particle, such as area or perimeter. The 
Quant imet usually ignores any holes within the measured area, but 
significant error is introduced when the perimeter is altered. 
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In considering the size of the speckle in the image, we find 
that the speckle arises from two causes [6] within our images. The 
first occurs in the recording of the hologram. Objects illumi- 
nated with diffuse coherent light will have a speckle pattern 
superimposed on them with a maximum speckle size d S p ec ^j e given 
by[6] 


X 

** speckle = (4) 

0 cos 0 

where X is the laser wavelength, is the angle subtended by the 
illuminating beam as seen from the object (see figure 6), and 6 
is the angle measured at the object between the axis to the 
center of the illuminated diffuser and the center of the image 
plane (i.e., the hologram). For the geometry of the present 
experiment, A=676 nm (a nominal wavelength between the recording 
laser and the reconstruction laser), 0 is approximately 62.4 
milliradians, and 6=180°, giving a speckle diameter on the holo- 
gram plane of 10.8 pm. This speckle will be recorded by the 
hologram and will be present in the reconstructed image. The 
primary means of reducing the speckle size is to increase the 
illumination aperture of the diffuser or to reduce the distance 
from the diffuser to the object. 

The second contributor to the speckle in the image is the 
imaging lenses. In the present system there are two imaging 
operations in the reconstruction process. The first is the forma- 
tion of the real image on the spinning mylar disk through the 
assisting lenses of the holocamera. The second is the imaging of 
the real image onto the face of the scanner tube through the 
microscope lenses. The maximun speckle diameter for a diffuse 
image that is imaged through a lens as depicted in figure 7 is [6] 


1.22XM 


‘speckle 


2 tan[ sin~l(NA)] 


(5) 


1.22XM 
2 tan[NA] 


( 6 ) 


where M is the magnification of the imaging system (given by the 
ratio of the image distance to the object distance) and NA is the 
numerical aperture of the lens. The numerical aperture of a 
microscope lens is approximately given by 

NA = (l/2f no ) (7) 

- tan -1 [d lens /2do] (8) 

where f no is the f-number of the lens* d 0 is the object distance 
from the lens, and dj eng is the diameter of the lens. Alternate 
forms of the equation for the maximum speckle diameter are 
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(9) 


1.22Xdj 

^speckle = 

d lens 


d speckle = 1.22X(l+M)f no (10) 

The first imaging operation is the formation of the real 
image on the mylar disk from the hologram through the assisting 
lens set. Although the lens set is involved in the imaging# the 
primary imaging element is the hologram itself. The hologram can 
be modeled as a lens with (approximately) unity magnification, 
the aperture of the lens is the width of the hologram as modified 
by the assisting lenses. Further work is necessary to complete 
the modelling of this imaging process. 

The second imaging operation is the imaging performed by the 
microscope objective onto the face of the scanner tube face. The 
maximum speckle size is given by equation 10 and is plotted vs. 
the f-number of the lens in figure 8 for a variety of magnifica- 
tions. Reduction of the speckle diameter to values below 3 or 4 
pm will require f-numbers below 0.4 for magnifications of lOx. 
For magnifications of lx, f-numbers below 4 are adequate. It is 
noted that typical f-numbers are in the range of 1.5 to 10 for 
this range of magnifications# indicating that magnifications of 
2x or so are consistent with reducing the speckle to values that 
are comparable to (or below) the value of the smallest expected 
particle. 

Generally the conclusion is that# to reduce the speckle 
size# large aperture optics are desirable. Once the speckle is 
smaller than the typical particle size# the logic test based on 
the feature size can be implemented in the Quantimet to extract 
particle information from the background. 

If the speckle cannot be made negligible compared to the 
particle size# other techniques must be used to reduce the con- 
trast level of the speckle relative to the particle. These tech- 
niques are based on the fact that the speckle pattern will shift 
if one or more of many variables are changed while the object 
position remains fixed. By averaging N images with differing 
speckle patterns, an reduction in the speckle contrast on the 
order of (N) -1 ^ is expected. With reduction of contrast the 
speckle level can be reduced below the grey level threshold set 
in the 1-D Auto detector module of the Quantimet. Variables that 
can be changed to cause a shift in the speckle pattern are: 

-Changing the random phase of the diffuse illuminating wave. 
This is done in our system by rotating the mylar disk that 
has the real image focussed on it. The speckle pattern 
changes during the integration time of the scanner tube. The 
improvement in contrast of the observed image (on the view- 
ing screen) formed with a moving disk over that of a sta- 
tionary disk is noticeable. The predicted improvement at 
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higher spin rates or longer integration times is less no- 
ticeable to the eye. 

-Changing the image aperture. Since the Quantimet requires 
equal strength images for accurate particle sizing, this 
technique would require a renormalization of the image amp- 
litudes. This might be possible when the system is in the 
advanced computer-controlled mode, but the technique holds 
little promise since one is limited by the limited dynamic 
range of the sensor to few images. 

-Changing the object position or, equivalently, changing the 
image screen position or changing the imaging lens posi- 
tion). This can be done within the depth of focus of the 
particle. For small diameter particles few independent 
images could be expected. 

-Computer superposition of images with the different diffuse 
illumination. This technique is similar to the rotating disk 
illuminator except that the superposition is done digitally 
rather than relying on the integration properties of the 
imaging tube. Limitations here are primarily in the memory 
capacity of the computer. Using the computer to store and 
manipulate the images also allows a variety of linear and 
nonlinear image processing algorithms to be used. 

Each of these techniques introduces complexity into the proces- 
sing of the image. The easiest to implement is the spinning 
diffuser in the image reconstruction. While the image is drama- 
tically improved when the disk is spinning compared to the image 
on a motionless disk, the image quality is not dramatically 
improved with increased rotation rate. The techniques other than 
computer processing the images appear to present insurmountable 
problems, leaving the computer processing as the most likely 
solution. Work has begun on implementing the computer-controlled 
configuration in an effort to explore this technique. 


SUMMARY 

The holographic technique has been successfully used to re- 
cord particles from strand burns and 2-D motor combustion. (Work 
will soon be done on expanding the hologram recording to a 3-D 
motor environment) The image particle count and sizing has suc- 
cessfully been done by hand but automatic techniques are prefer- 
red. The diffuse illumination required in the recording process 
to avoid phase fringes due to thermal effects has led to the pre- 
sence of speckle in the images. Efforts have been made to analyze 
the size of the speckle to reduce the maximum speckle diameter. 
The remaining speckle is reduced by using a spinning diffuser in 
the reconstruction process. The computer-processed images appear 
to have the most liklihood of further reduction f the speckle 
contrast. This method is currently under investigation. 
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Figure 1. Schematic of 2-D Motor 
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Figure 2. Schematic of Reconstruction 
Viewing Method (after Ref. 3) 



ORIGINAL PAGE 53 

OF POOR QUALITY 



(b) WGS-ZrC 


Figure 3 Photographs of Typical Reconstructed Images 
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Figure 4 Block Diagram of Qoantimet 720 
EnhancedBaeic Configuration 
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Figure 5 Block Diagram of Semiautomatic 
Image Reconstruction Setup 
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Figure 6. Geometry for Calculating Maximum 
Speckle Diameter Due to Diffuse Illuminator 
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Figure 7. Geometry for Cal oule ting Maximum 
Speckle Diameter in Imaging 
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SPECKLE SIZE, micrometers 



Pigure 8. Maximum Speckle Diameter vs. 
Leas f-auaber aad Magaif icatioa 
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